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Cosmic Ray Propagation 
C. E. FICHTEL and D. V. REAMES 
Goddard Space Flight Center 
Greenbelt, Maryland 

ABSTRACT 

The general problem of the propagation of cosmic rays from their source 
until they are observed in the vicinity of the earth is reviewed. The effect of 
matter on the cosmic ray composition is considered in detail and the effect of 
diffusion is treated in the depth necessary to study some models of the cosmic 
radiation not previously considered in this manner. The results of these calcu- 
lations show that there are some aspects of the experimental observations, such 
as the light to medium nuclei ratio as a function of energy, which are not con- 
si stent with any equilibrium model of the cosmic radiation which does not include 
a rigidity dependence in the mean free path. This conclusion is independent of 
the source spectrum assumed and is valid for a large class of source distribu- 
tions. A series of possible explanations for the discrepancy between the ob- 
served experimental data and the theoretical predictions of the equilibrium pic- 
ture are considered, and shown to be unlikely. Faced with this dilemma, 
attention is turned to a simple non-equilibrium model, consisting of one fairly 
close single source superimposed on a general background. This picture was 
chosen since it seemed the simplest first step from equilibrium. The calcula- 
tions show that this picture leads to considerably better agreement with the ex- 
perimental data, especially the light to medium nuclei ratio as a function of 
energy/nucleon. At least one irreconcilable piece of experimental data remains, 
namely the one measurement of the Fluorine to Oxygen ratio. 
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I. INTRODUCTION 

The origin and history of the energetic cosmic rays which continually bom- 
bard the earth is a subject of considerable interest because of its relation to 
so many fundamental aspects of the galaxy and possibly the universe, including 
the energy content, the interrelationship of matter, magnetic fields, and cosmic 
rays, and the nature of the objects which are able to accelerate particles to such 
high energies in very great numbers. In attempting to obtain a fuller understand- 
ing of these questions, one of the fundamental problems which must be examined 
is the propagation of the cosmic rays from their source to the vicinity of the 
earth where they are observed. 

It is generally believed that after leaving their sources, cosmic rays diffuse 
through interstellar space with their motion being controlled and made random 
on a large scale by the magnetic fields. In their passage through interstellar ma- 
terial at least two mechanisms affect the cosmic-ray composition and energy spectra; 
these are fragmentation produced in nucelar reactions with the interstellar 
material and Coulomb interactions including ionization energy loss. Additional 
processes which might affect the energy dependence of the relative composition 
include Fermi acceleration in collisions of the cosmic-ray nuclei with magnetic 
irregularities in clouds and a rigidity dependent escape from the galaxy, or stor- 
age region, if the cosmic rays are limited to some region and do not pervade all 


of the universe. It is also possible that cosmic rays have passed through some 
material before leaving the source region, where the above processes may also 
occur. 

In an earlier paper (Fichtel and Reames, 1966)! the theory of energy de- 
pendent propagation of cosmic rays through interstellar space was developed 
under the following assumptions: (a) The source energy/nucleon spectra of all 
multiply charged nuclei have the same shape, at least above 1 00 MeV/nucleon. 
(Note that, since all of the multiply charged nuclei that will be of interest at the 
source have nearly the same charge to mass ratio, this effectively permits the 
spectra to be both velocity and rigidity dependent.) (b) The relative abundance 
of He 3 and light nuclei (3 < Z ^ 5) at the source are negligible compared to He 4 
and medium nuclei respectively, (c) The average interstellar potential path 
length from the source to the earth is independent of the energy /nucleon of the 
particle, (d) The potential path length distribution was reasonably smooth and 
did not contain a high percentage of very long or very short paths. In this same 
article, alternate approaches and models suggested in the literature were also 
examined. 

Subsequently, new measurements on the fragmentation cross sections for 
cosmic rays interacting with interstellar matter have been made, and, in some 
instances, their values are markedly different from those previously assumed. 
In light of this and new measurements on the relative abundances of cosmic 
rays, the calculations were repeated, and the results presented in a recent paper 
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(Reames and Fichtel, 1967) together with a consideration of the problem of pos- 
sible interstellar acceleration. The relative abundances deduced in that paper 
indicated that there seemed to be no simple way of explaining the experimental 
results particularly the relative abundance of the light (3 ^ Z ^ 5) and medium 
(6 $ Z< 9) nuclei within the limits of the assumptions. 

The application of the same propagation method to different types of media 
including one which more nearly simulated the source region was accomplished 

3 

by Durgaprasad. He considered an ionized medium, one consisting of a compo- 
sition which approximated that of a supernovae as nearly as possible, and a 
partially ionized region with a small percentage of helium as well as hydrogen 
to better approximate interstellar space. The results of his work showed that 
within the range of interest the shape of the rate of energy loss curve was 
virtually the same in all cases, but that the magnitude of the rate of energy loss 
was greater for ionized regions. Also, there was no significant difference be- 
tween the rate of energy loss in pure hydrogen and the combination approximating 
interstellar space. Thus, in this work there is no need to be concerned about the 
exact interstellar composition. If the cosmic rays have spent a significant per- 
centage of the material path length in an ionized source region there is an ef- 
fect, which can generally be summarized as the suppression of heavier elements 
at lower energies. 

One approach which goes beyond the assumption of a single path length for 
all particles involves assuming that the cosmic rays are in equilibrium at a given 
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point and further that there are no spatial variations or second order energy 

2 

effects. It has been shown previously that even with the inclusion of a Fermi ac- 
celeration effect this approach will not give satisfactory agreement with the 

4 

experimental data. Cowsik et al. (1966) nave suggested adding a term to the 
equation for loss of particles atapoint which is intended to represent the escape 
from the galaxy in an equilibrium model. Actually, an escape of this kind should 
be included as a boundary condition in the solution of the complete equation and 
not as an additional term. Once boundary conditions must be applied, an analy- 
sis involving solving a differential equation exactly or one such as that to be 
presented in this paper is appropriate. The mathematical results derived from 

introducing a term for loss at a point are given by Cowsik et al.? and in Reames 
2 

and Fichtel . In the latter paper the predictions are shown to disagree with the 
experimentally observed light to medium nuclei ratio. 

In comparing experimental results to theoretical predictions, the problem 
could be complicated by the fact that the local solar modulation is not yet known, 
but the general belief is that it probably depends only on the velocity and 
charge-to-mass ratio of the particle. Therefore, although nuclei of the same 
charge-to-mass ratio, but different charges will lose energy at different rates 
in interstellar space, the fluxes of these particles will be modulated in the 
same way, thereby permitting the separation of modulation effects from inter- 
stellar energy loss and fragmentation effects. 


4 
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In this paper, we wish to review and extend the previous work, first by 

summarizing the fundamental considerations and expanding parts of the previous 

work to show more clearly the nature of the problem, and secondly by pursuing 

the implications of the results and the possible reasons suggested in the 
2 

previous paper for the deviation of the experimental results from those ex- 
pected on the basis of the earlier calculations. Other work on the same subject 
will also be reviewed at the appropriate points in the discussion. The con- 
clusions reached after an examination of various alternatives will be presented. 
A two source model wherein one recent local source is superimposed on a gen- 
eral cosmic ray flux will be given special consideration. 

■* 

II. PROPAGATION THEORY 

This section will be divided into three parts. The first part will be devoted 
to a brief summary of the problem of the effect of the travel of particles through 
matter and the method used to determine the relative abundances and energy 
spectra after travel through any amount of matter. The second part will involve 
a difecussion of the diffusion problem, which must be treated more explicitly now 
in ahticipation of the need to consider more complex cases than have been ex- 
amined formerly. 

A. Effect of Matter 

The present approach to energy — dependent propagation of cosmic rays 
through interstellar space is based on the one described in detail previously , 1 
which includes both energy loss and fragmentation. The fundamental transport 
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equation used there is : 

^ [»!<*,. X)j, (£,.»)] * »i(E,.X)j i rE,,x)/A, 

= «i(E.,X) £ j fc CE..*VA tl (E,.X) (D 

k>i 

where j i (E s , x) is the flux per unit energy /nucleon of i-type particles of energy/ 
nucleon E after propagation through x g/cm 2 of material given their initial 
energy/nucleon to be E s , w i (E s , x) = (dE/dx); for these particles, A. is the loss 
mean free path, \ i is the mean free path for production of i-type particles from 
k-type particles, and E s is the energy/nucleon at the source. The specific assump- 
tions mentioned in the introduction were discussed in the earlier article 1 in detail. 
Since the reasons for believing the first two have not changed, these assumptions 
will be kept; however, the third assumption will be examined later, and the fourth 
in sections IIB and III. 

The calculations are made in steps of 0.02 g/cm 2 and the individual ele- 
ments from Helium through Oxygen and the charge groups 9 ^ Z ^ 19 and Z > 20 
are considered. The dominant cross sections used in the fragmentation process 
are shown in Table I. Most of the cross-sections listed are taken from experi- 
mental measurements 5 “ 9 . Note that the low loss cross section at low energies 
for VH nuclei reflects the fact that most of the total reaction cross section re- 
sults in the production of nuclides of only slightly lower change; the probability 
of producing a nucleus 4 with 3 < Z<20 for a nucleus with Z = 26 is extremely small 
in this region. Loss cross sections for H nuclei are similarly effected. Since 
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oxygen represents the heaviest single element considered, its loss cross section 
is more nearly the measured total interaction cross section as in the case for 
lighter elements. These cross sections depend approximately upon A 2/3 . The 
cross sections for production of light nuclei shown in Table I are now reasonably 

5-8 

well known. 

B. Diffusion, Source Distribution, and Boundary Conditions 

Previously Fichtel and Reames 1 noted that because of the smooth varia- 
tion in the flux ratios from 1 to 6 g/cm 2 , the assumption that all cosmic rays 
have traversed a given amount of material, x o » will give essentially equivalent 
results to assuming a reasonable path length distribution which has an average 
value of X Q . That this statement is basically correct is evident also from the results 
of Balasubramahyan, et al.; 10 however, we now wish to consider at least some cases 
where there is a substantial probability for very long and, more importantly, 
very short path lengths. Therefore, it is appropriate at this point to review the 
problem of the diffusion of cosmic rays. 

The flux Fj of particles of energy between E 0 and E 0 + AE at the observing 
point, r 0 , at r is given by the expression 

= f <*>'. f <l,.r..T) (2) 

v T j > i 

s 

where the S . are the actual sources intensities, 
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Table I 

Dominant Fragmentation Cross Section (mb) vs Energy (MeV/nucleon) 


E 

30 

50 

70 

100 

150 

200 

300 

500 

800 

1000 

4000 

10,000 

Reaction 
VH loss 

* 

0 

0 

0 

0 

0 

0 

0 

7 

26 

44 

220 

400 

H loss 

36 

67 

88 

105 

117 

124 

133 

143 

153 

160 

200 

300 

0 loss 

570 

450 

300 

290 

290 

290 

290 

290 

290 

290 

290 

290 

£ 

T 

o 

138 

171 

154 

125 

98 

90 

84 

84 

86 

89 

90 

90 

o - c 

136 

96 

75 

61 

55 

52 

51 

50 

50 

50 

50 

50 

VH - B 

0 

0 

0 

0 

0 

0 

0 

2 

4 

6 

24 

32 

H -> B 

0 

0 

3 

7 

11 

13 

16 

19 

22 

24 

24 

24 

O - B 

3 

25 

38 

32 

25 

25 

25 

25 

25 

25 

25 

25 

C - B 

114 

128 

118 

98 

73 

63 

54 

50 

49 

49 

49 

49 

O - Be 

0 

2 

5 

7 

8 

8 

9 

9 

10 

11 

11 

11 

C - Be 

5 

18 

14 

11 

10 

9 

9 

11 

12 

12 

12 

12 

VH -* Li 

0 

0 

0 

0 

0 

0 

0 

1 

3 

5 

15 

22 

H - Li 

0 

0 

0 

2 

6 

8 

12 

18 

23 

25 

30 

30 

O - Li 

0 

2 

8 

16 

22 

26 

30 

35 

36 

36 

36 

36 

C - Li 

4 

18 

20 

19 

16 

14 

12 

12 

12 

12 

12 

12 

B - Li 

36 

35 

32 

30 

30 

30 

30 

30 

30 

30 

30 

30 

Li loss 

250 

170 

120 

110 

110 

110 

110 

110 

110 

110 

110 

110 


*A zero implies only that the cross section is less than 1 mb and hence negligible for the calcu 
lation. 


is the propagation function for the i'th species of energy E 0 at r 0 and t arising 
from energy E s at the source r s at time r before t. if the medium in which the 
particles are diffusing is isotropic with a constant density 8 , or can at least be 
approximated by an average density, the quantity t - r is given by the equation 
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(t - r) = 1/ 'c/3 = x/S c/3 (3) 

where -£ is the path length in conventional units and x is the path length in 
density • length (e.g. g/cm 2 ). Further, if S j is independent of time, or at 
least on the average over the time scales of interest, then eq. (2) becomes 

dl 

j>i 

In the last equation, Y. i includes both the effects of fragmentation and energy 
loss, discussed in the previous section, and the diffusion function. These two 
effects are strictly separable if the diffusion function is not a function of the 
type of particle or its rigidity. In fact, it is still a fair approximation to con- 
sider them as separable as long as there is not an appreciable change in the 
diffusion characteristics over the range of rigidities involved in the transit of a 
given particle from the source to the earth. Under the assumption mentioned, 

Y 3 1 is the product of a diffusion propagation term D {l , E 0 , r s , r Q ) and a material 
propagation term M.j ( E s , E o>^)- 

Equation (4) then becomes 

F,(E„.r 0 ) = dtD(-C,E 0 ,r > ,r o )£]»l jl (E s ,E 0 ,-t)S j <B s ,r,) (5) 

V s 0 j>i 

The last term in eq. (5) represents the result of propagation through % g/cm 2 
of material and is, therefore, the solution of eq. (1) described in part HA, i.e. 


Fi(E 0 . r o)=J d3r s J 


L 


V( E s- r s - E 0 -r 0 ^) Sj (E s ,r s ). 


(4) 
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i i ( E o-^) M ji( E s- E o^) Sj(E 8 ,r B ) (6) 

j — i 

Eq. (5) may now be written as 

j.( E 0 ,^)d^ (7) 

D is determined in the manner to be described now, and the integral over d-t is 
accomplished by a means of a computer in steps of 0.02 g/cm 2 . 

Eq. (7), in effect, simply states that theenergy spectrum deduced for each type 
of particle for a given path length i is multiplied by the relative probability of 
occurence of that path length, and then this product is integrated over all pos- 
sible values of the path length. 

In order to obtain an explicit expression for D d 3 r, which is in effect a po- 
tential path length distribution for cosmic rays, specific models for the origin 
of cosmic rays must be considered. At one extreme is a single source region, 
perhaps the galactic center, and at the other perhaps is the concept of cosmic 
rays spread uniformly throughout the entire galaxy or even the universe and 
fed by sources spread throughout this whole region. An intermediate picture 11 , 
which is often thought to be closer to the real situation, is one in which the 
cosmic rays observed in our galaxy are the result of many supernovae which 
are randomly spaced in time about a hundred years apart and probably concen- 
trated towards the central part of the galaxy and in the spiral arms. In the first 


^(Eo-r,) 
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part of the work that follows we shall assume that in the past, or at least that 
part in which a large part of the cosmic rays observed are produced, the proba- 
bility of cosmic ray production was independent of time, because of the time 
scales involved. This assumption implies that all of the situations which we are 
describing are equilibrium ones, which is reasonable for the cosmic ray case. 

An alternative concept is that they were all produced at one point in the past 
presumably the origin of the galaxy. The latter theory has the difficulty that the 
density of matter is such that the composition of the cosmic rays would probably 
be very different from that observed due to the very long path length. 

Another possibility is that a single relatively close supernovae produced 
a significant fraction of the cosmic rays. This assumption will be explored in 
part IV, and the results must be reconciled with the considerable evidence sug- 
gesting the constancy, to within a factor of two of cosmic rays over time periods 
which are long compared to the time needed to traverse the average path length. 

A few mathematical models related to the problem under consideration will 
now be examined and general conclusions will then be drawn from the results 
of the calculations. In looking at the problem of diffusion in the galaxy, we are 
faced with the difficulty of having little idea of an appropriate model. Therefore, 
we shall examine several extreme cases to try to place constraints on the type 
of processes which might be occurring. 

The consideration will begin with the diffusion from a point source into 
a uniform medium of finite extent. This will be followed by the cases where 
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the medium has an infinite extent and where there are a distribution of sources 
The distribution of potential path lengths from a point source to an observation 
point a distance r away can be obtained readily from diffusion theory , and is 
given by the expression 


JL±T " sin exp AlAlM') 

2 d 3 r rrl ' d • \ 3d 2 / 


where d is the distance to the edge of the medium, \ is the mean free path be- 
tween scatters, -£ is the path length, and N is a normalizing constant which de- 
pends on the intensity of the source. To obtain the expression in eq. (8), the flux 
at the boundary is actually assumed to be zero. When this approximation is not 
made a much more complex form is obtained; however, for the purposes of this 
work, the expression above is quite adequate. This formulation also assumes 
that the source has been uniform in intensity throughout the past. The calcula- 
tion which follows immediately now can therefore be applied to the cosmic ray 
case if the time between outbursts in a region to be considered is such that vAt 
(= A -£) is small compared to the average distance. Taking the particular ex- 
ample of supernovae in our galaxy At - 100 years for the galaxy as a whole, or 
from 10 to 10 times this value for a smaller region. If v = c, and there is 1 
atom/cm 3 , A^= 2 x 10 -2 to 2 x 10 -3 g/cm 2 compared to an estimated average 
path length for cosmic rays of the order of 4 g/cm 2 . The value of -i depends 
somewhat on the specific model, but this value is quite adequate for this present 
discussion. 
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Notice first that as ■i becomes sufficiently large so that ( 77 2 /3d 2 ) » 1 

eq. (8) becomes: 



N 

2d 3 



(9) 


Any distribution considered then in a bounded medium for sufficiently large 
will ultimately approach an exponential character. For the cosmic ray case 
the condition (7r 2 \-£/3d 2 ) >>1 for most of the range of i values of interest would 
imp ly the following: From the experimental data on the light to medium ratio we have 
already mentioned that -t is about 4 g/cm 2 or 2.5 x 10 24 cm under the as- 
sumption of 1 atom/cc. The radius of the galaxy is about 4 x 10 22 cm, which 
is the largest trapping region that seems reasonable to consider. Thus, k must 
be much larger than 2 x 10 cm if eq. (8) were to be a reasonable approxima- 
tion over most values of l . Since this would imply X. was only one or two orders 
of magnitude smaller than the galaxy and of the order of the thickness of the 
galactic plane, it seems unlikely that the equation (7T 2 k'i/3d 2 )» 1 holds for -i 
values of interest. 

Returning to eq. (8) if d » r and d 2 »/^ , the sum can be replaced by an 
integral. Upon completing the integral the following equation is obtained: 


DA't” 


d-*co 


3 3/2 N 

47t 3/2 K 3/ H 3/2 


exp 


A3rA 

\iu) 




(10) 
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Since this equation is increasingly close to being exact as d -co , it is the expres- 
sion for an infinite isotropic medium. 

The second condition probably holds for the galaxy for l values of interest 
and the first will hold for sources which are not too distant. If the first condi- 
tion does not hold very well, the physical result is not expected to be too different. 
Notice also that as r becomes large Dj decreases very quickly as the square 
of r 2 in the exponential, except for values of K, much greater than those of 
interest. As will be seen except for a very unusual source distribution the 
case of large r, where r - d , will not be of interest. 


Consider now the case where eq.(10) is valid, but instead of one source 
there is a distribution of sources, with a density given by p, where p includes 
both the effect of varying source strength and source density. Here, we are 
assuming that there is no source which has been so close to the observer or so 
strong as to cause it to stand out over the others. If this assumption is true, 
the sources can adequately be represented by a distribution function p . If p 
is normalized such that JJJ pdv = l, then the distribution of path lengths will 
be given by the expression 


JDM. d 3 r s 



3 3/2 N A l 

877 3 / 2 \ 3 / 2 4 3/2 


p exp 



r 2 s i n0d0d<I> 


( 11 ) 


Letting 


z 


t/v 4X^/3, 


( 12 ) 
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yields 


J DA't d 3 r s 



joz 2 exp (-z 2 )dz sin^d^dtf*. 


(13) 


Let us now consider certain special cases for the function p, assuming for the 
present that k is not a function of position. 


Case (a): Suppose the sources are uniformly distributed over all space 
with p = p 0 . Then (12) becomes: 


NA-tp f° 

jDd 3 r s A'fc= 47 t I z 2 exp (-z 2 ) dz 

s 7 7- 3 /2 Jq 

(14) 

jDd 3 r A4 = (N / o 0 )A^ . 


Notice that this result does not depend on -i and that it is one example 
wherein there are no spatial or energy gradients of second order so that when 
it is inserted in eq. (7) the result is identical to the approximate equilibrium at 

2 4 

a point solution discussed earlier ’ . 

Case (b): Suppose next that the sources are limited to some finite volume, 
and, to obtain an explicit expression, assume that p- p 0 for r < r a and p = 0 
for r > r a . Then eq. (12) becomes 

3 r s A'E = q 477 f z 2 exp(-z 2 ) dz ^5 

77 3/2 4) 
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where 

a = r f / 4\-t/3 . 

This last integral can be integrated by parts to obtain 

jDd 3 r s J dzj . 

The function given in eq. (15) is plotted in fig. 1, using the tables given 
by Jahnke and Emde 13 to evaluate the integral. 

There are now four specific distributions that have been developed, and the 
implication of these results will be examined since most other possibilities 
will lie in between these cases. For example if the distribution of sources in 
the medium is not centered about the observing point, the potential path length 
distribution will appear somewhat like a superposition eqs. (10) and (15). If 
the sources all lie appreciably away from the observing point the potential path 
length distribution will appear similar to eq. (10) only flatter and broader. 

Thus far, k has not been investigated in detail. In principal k can be a 
function of both the position in space and the rigidity of the particle. In fine 
structure A will certainly vary strongly with position; whether or not it does on 
a larger scale is not an easy question to answer, except that it probably in- 
creases toward the outer boundary of the galaxy. Eq. (8) shows that this effect 
would increase the tendency towards an exponential distribution in ■i at a smaller 
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l value. Also, if it is assumed k is equal to A- 0 (r/r 0 )^ eq. (10) is replaced 

by 


r 


exp 


3r^r 2 " 

(2-/3) 2 \/J 


( 16 ) 


Further considerations of this type will not be pursued because as we shall see 
later, the resulting predictions for cosmic ray spectra and abundances are very 
insensitive to changes in forms of this type. 

A much more significant consideration is the possible variation of k with 

rigidity. Clearly it is unlikely that all cosmic ray particles will be affected the 

same way by the magnetic irregularities in space. In general, a distribution 

in scale sizes of the irregularities would be expected, and, since the particles 

of lower rigidity would then be expected to encounter more significant deflec- 

1 2 

tions, their mean free path would be longer. However, as Parker has pointed 
out in relation to solar particles, it is also possible for low rigidity particles to 
follow field lines along a large kink which would affect a high rigidity particle. 
Since essentially nothing is known about the dependence of mean free path on 
rigidity except that once the rigidity is large enough, the path length will de- 
crease with energy, a trial function will be used with the aim of seeing 
what effect it has. Assume that the mean free path was given by the ex- 
pression k = \ 0 + AR, which states that k approaches a constant at low rigidities 
and then increases in proportion to the rigidity at large values. The appli- 
cations of the mathematical expression for k is restricted by the limitation 
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mentioned earlier in the section that D cannot depend strongly on E Q , or it is 
not valid to separate it from the material propagation term. (See the discussion 
following eq. (4).) 

Notice, however, that the eq. (14) for sources spread uniformly through- 
out space does not depend on \ and hence cannot depend on rigidity. 

HI. EQUIUBRIUM MODELS 

A. Comparison to Experimental Results 

In this section, the resulting relative abundances of the cosmic ray nuclei 
will be studied as a function of the source spectral shape and the potential path 
length distributions developed in section II for cases of equilibrium. This study 
will be supplemented by a discussion of other equilibrium models suggested in the 
literature. 

As in the earlier paper, 1 we shall choose three trial source energy spectra 
ranging in steepness from a power law in total energy to a power law in kinetic 
energy. The specific expressions used are 

dj/dw = C a W " 2 - 5 ( a) 

dJ/dR = C b R- 2 - 5 (b) 

dJ/dE = C c E" 2 - 5 (c) 

where W is the total energy, R is the rigidity, and E is the kinetic energy per 
nucleon. 
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For comparing the results here, the light to medium nuclei, the helium to 
medium nuclei and the heavy to medium nuclei ratios were chosen to be dis- 
played because the first shows some of the more significant effects since the 
light nuclei are secondaries and since the latter two ratios are typical of many 
of the others. 

For comparison, the results for the case where the potential path length 
was a delta function, i.e. the common fixed path length approximation, are shown 
first in figs. 2,3, and 4. Here the path length was adjusted to give a high-energy 
L/M value consistent with the data. The path length distribution will in general 
be normalized in this way whenever there is a free parameter. Notice first that 
the better cross-section data now available no longer indicate that any increase 
in the light to medium nuclei ratio should be seen at energies in the neighbor- 
hood of 200 to 500 MeV. Therefore, if the increase seen by some experimenters 
is real, it is not a fragmentation effect. This effect does not depend on the 
model as we shall see. Notice also that L/M is essentially independent of 
energy and that He/M rises very rapidly at low energies for steep spectra. 

Turning now to the case of a point source expanding into an infinite isotopic 
medium, eq. (10) applies, and the results of the calculation outlined in the 
previous section lead to figs. 5,6, and 7. These graphs show that there is little 
difference from the previous case as expected since eq. (10) states that there is 
not a heavy weighting of very short nor very long path length. 
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The next case considered in section IIB was the uniform distribution of 
sources over an infinite medium, which led to a potential path length independent 
of i given by eq. (12a). Here, unlike the other cases, there is no adjustable 
parameter which can be used to normalize the result to the high energy L/M 
value. Hence, since the calculated light to medium nuclei value shown in fig. 8 
disagrees with the observed one at high energies, this case can be rejected at 
once, when k is a constant. 

An important more general statement can now also be made. It was pointed 
out in section IIB, that the same approach to a solution, namely simply inte- 
grating over all path lengths with equal weight, also applied to any equilibrium 
situation wherein there are not serious gradients in space or energy. Hence, 
the experimental results on the light to medium nuclei ratio, also exclude thi£ 
class of model. One important set in this class is any universal model for the 
cosmic rays except those which somehow would predict substantial spatial or 
energy gradients. 

Eq. (15) referred to sources limited to a region around the observing 
point. Figs. 9 and 10 shows that the light to medium nuclei ratio decreases 
markedly at low energies whereas the helium to medium nuclei ratio remains 
fairly constant. These characteristics appear in cases where there is a 
reasonably high probability for short path lengths. The reason is that the low 
energy spectra changes rapidly with distance traversed because the rate of 
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energy loss is an increasing function of decreasing energy. The result is that 
the steep initial spectrum is bent downward quickly; so, if there is a relatively 
high probability for short path lengths, the low energy region reflects the source 
characteristics quite strongly. Thus, a relative low abundance of light nuclei is 
seen, and the helium to medium nuclei ratio has not had a chance to increase 
much as it will with additional passage through matter because of the greater 
rate of energy loss of heavier nuclei. 

Another case of possible interest is the one in which all of the sources 
are at a distance which is greater than some value. This possibility can be 
approximated by assuming that sources are uniformly distributed beyond some 
radius r from the observing point. The resulting path length distribution 
would then be obtained by subtracting some multiple of eq. (15) from eq. (14). 
Regardless of the assumptions about the value of \, except that it is a con- 
stant, the L/M value will now be even larger at high energies than for the case 
of uniformly distributed sources since the average path length is even larger. 
Thus, this model by itself is not in agreement with the experimental result. 

The results of the above show that none of the above distribution give a very 
satisfactory agreement with experimental results on the light and medium nuclei 
unless the increase at low energies due to the modulation effect were much 
stronger than expected. 

Considering next the case where \ is possibly a function of rigidity. Fig. 

11 shows the result for a point source expanding into an infinite isotopic medium 
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with X = a (1 +bR ), i.e. essentially a constant at low rigidity and proportional to 
rigidity for large values. The light to medium nuclei ratio variation is in reason- 
able agreement with experiments over the measured range as are the other 
ratios, but it should be remembered that this form would predict that the light 
to medium nuclei ratio would decrease at high rigidities. In general, it seems 
possible that X should ultimately increase with rigidity, and, hence, if the cos- 
mic ray sources are in the galaxy, then the light to medium ration should de- 
crease at very high energies. However, the variation with energy could be 
much smaller than that predicted by a X dependence as extreme as the one 
chosen above to obtain agreement at low energies. Further, such an extreme 
rigidity dependence implies a very steep source spectrum, j s ~ R _1 5 
instead of j s ~ R" 2 ' 5 . Using an approach wherein it is assumed that the vari- 
ation of the light to medium nuclei with energy is in fact due to a variation of the 
mean potential path length with rigidity, Apparao, 14, 15 and Biswas, et al. 16 ’ 17 de- 
duced a variation of path length with rigidity which peaks around 200 to 500 MeV/ 
nucleon, similar to the result here. Kaplon and Skadron 11 suggested that this 
dependence of the path length on energy /nucleon may be the result of a rigidity 
dependent mechanism at the source. It could also be due to a rigidity dependent 
path in space as mentioned earlier. It seems difficult to explain why the rigidity 
dependence needed to obtain agreement with the experimentally observed light to 
medium nuclei ratio should exist. Specifically why should there be a strong rigidity 
dependence between 0.3 and 1.5 BeV/nucleon and apparently very little above, 
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assuming the light to medium nuclei ratio is approximately constant above 1.5 
BeV/nucleon. If the light to medium nuclei ratio should continue to decrease 
markedly with increasing energy/nucleon, then a very flat source spectrum is 
implied as mentioned before. 

Another set of components of interest are deuterium and He 3 . These 
elements are presumably nearly absent in the source; so whatever fluxes are 
seen are presumably secondary. Experimental data 18-22 on these nuclear 
species only exist in the low energy (< 300 MeV/nucleon) region. The 
interpretation of the He 3 data is complicated by the fact that the He 3 charge to 
mass ratio is appreciably different from that of the parent nuclei, and, hence, 
He 3 is modulated differently from He 4 in the solar system in a way which is 
not yet well determined. Ramaty and Ligenfelter 23 have considered the problem 
in detail, actually using the combined data to estimate the solar modulation, 
and found the data to be consistent with an average interstellar path length of 
4 ± 1 g/cm 2 . Since their work is complete in itself and agrees with that of most 
others ' who have considered this particular problem, the calculation re- 
lated to these secondary particles will not be repeated here. Rather it will 
simply be pointed out that this result is consistant with the results presented 
already for the heavier nuclei in the low energy region except for spectra which 
are quite steep, and, as in the case of the heavier nuclei, the result is generally 
not sensitive to the particular diffusion model chosen. 
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A serious problem arises, however, in relation to the one measurement 
made so far on the relative abundance of fluorine. After the cosmic rays have 
passed through 3 to 5 g/cm 2 of material, the fluorine abundance should be 
about 1 or 2% of the abundance of oxygen or carbon due to fragmentation of 
elements heavier than fluorine. The appropriate cross section are not really 
known, but by analogy to similar reactions such as F 19 (p, pn) F 18 a reason- 
able estimate can be made. The principal cross section is the one for the 
production of fluorine from neon, which is about 40 millibanks at 100 MeV/ 

25 

nucleon and decreases towards higher energies. Comstock et al. found 
experimentally a fluorine to oxygen ratio of £ 1.4 x 10 " 3 for an energy interval 
around 10 2 MeV/nucleon. Comstock et al. 25 and Burbidge et al. 26 both suggested 
a two component model to explain this result, one component of the cosmic rays 
having gone through 0 g/cm 2 and the other through 20 g/cm 2 . These two com- 
ponents are then added in such a way to give the correct light to medium ratio. 
Using the procedure outline in section IIA, the composition at 20 g/cm 2 was 
calculated and added to the composition at 0 g/cm 2 in proportions which gave 
the light to medium ratio observed experimentally at 100 MeV. For an initial 
spectra wherein J ~ R 1 5 a fluorine to oxygen ratio of about .01 was obtained. 
Similar results were obtained for J ~ W -1 ' 5 and J E~ 1 ' 5 . Further, a lithium 
to boron ratio of about 5:3 is obtained compared to the experimentally measured 
value of about 2:3. Finally, if the two spectra are added to give the correct L/M 
ratio at low energies, a result inconsistent with the data is obtained at high 
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energies because the 20 g/cm 2 spectrum is much flatter and hence strongly 
dominates at high energies. Thus, when examined in detail this model also 
seems not to lead to predictions consistent with experimental results. 

Various other assumptions including those already discussed here will 
also not lead to a fluorine to oxygen ratio which is nearly as low as the experi- 
mental result of Comstock et al. Whether the fragmentation occurred in the 
source or in the interstellar medium is also irrelevant. There are still two 
possible ways in which the result could be consistent with the fragmentation 
concept of the formation of light nuclei. Firstly, the assumed cross sections 
could be wrong; secondly, the experimented measurement of the upper limit of 
fluorine could be wrong. Therefore, it would be extremely desirable to measure 
the neon to fluorine cross section as well as the production of fluorine by heavier 
elements. Also in view of the very unique character of the fluorine abundance 
measurement, confirmation by a second experiment would be desirable. 

B. Discussion of the Apparent Failure of the Equilibrium Model 

The failure to achieve agreement between the prediction of the models described 
thus far in this work and experimental results may possibly result from several 
features. In the following paragraphs we shall examine what to us appear to be 
the most obvious potential difficulties with the present approach and hence pos- 
sibilities for the discrepancies, but first it is worth reemphasizing that the 
alternatives are being considered below because there no longer seems to be 
any way to resolve the dilemma within the framework already discussed. 
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(A) Propagation through material in interstellar space or in the source may 
not be the dominant process for producing light nuclei and/or they maybe emitted 
directly by the source. This possibility has been considered from time to time , 27,1 
but, to our knowledge, no formal treatment in the literature exists. There is no 

3 

means known to us of forming the nuclei such as Li, Be, B, He , F during the 
acceleration process except by fragmentation. Also, it is not yet clear how the 
rather unique variation of the light to medium nuclei ratio with energy observed 
experimentally could be derived theoretically. 

(B) There may be further inadequacies in the cross section data. This 
possibility, although once of great concern, no longer seems as significant, due 
to the extensive new data on most of the significant cross sections. There are 
still some important gaps in the data, such as the production of fluorine. 

(C) The source spectra for particles of the same charge to mass ratio 
may not be the same. The theoretical reasons for believing similar source 
spectra exist at the source were outlined previously . 1 The existing experi- 
mental evidence for helium, medium, and heavy nuclei is roughly in agreement 
with the assumption of similar source spectra and an average interstellar path 
of a few g/cm 2 . This last statement is certainly true of carbon and oxygen the 
principal parents of the light nuclei. Finally, there is no known way, consistent 
with the observed experimental data, for the spectra to differ in a way which 
would significantly help in resolving the dilemma of the light to medium nuclei 
ratio. 
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(D) The observed features of the galaxy are considerably more complicated 
than those of the models . A reexamination of the models show that additional 
complications are not likely to resolve the problem unless they either (1) intro- 
duce a very unlikely rigidity dependence for the mean path length for scattering 
or (2) introduce a very unique rigidity dependent trapping at the source. 

(E) The energy spectral measurements on the light or medium nuclei and other 
species are wrong. Although this may be a possibility due to the difficulty of charge 
and energy measurements in the intermediate energy range, for the purpose 

of this work, we shall assume they are correct until subsequent work indicates 
otherwise. 

(F) Solar Modulation: There is almost certainly a rigidity dependence in 
the solar modulation, but its magnitude and form are not yet certain. Since light 
nuclei have a slightly higher rigidity for a given velocity, the effect of the rigidity 
dependent modulation is generally to increase slightly the light to medium ratio. 

23 

Estimates of the degree of modulation by Ramaty and Lingenfelter , Gloeckler and 

28 2 9 

Jokipii , and Durgaprasad et. al. seem to indicate that the increase is at most 

a 2 to 6% effect around 300 to 500 MeV/nucleon which is not enough to affect the 
previous discussion significantly, and would, for the most part, be hidden in the 
experimental errors. Under extreme assumptions, e.g. assuming the modulation 
function as of the form exp (-C//^R) down to low energies, the modulation effect 
could increase the light to medium ratio by 5 to 15% at about 100 MeV/nucleon, 
but more likely the rigidity dependence of the modulation decreases at low ener- 
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gies making the effect quite small. If, in the extreme, the modulation were that 
necessary to produce the low energy cosmic ray spectrum resulting from the 
supernovae theory proposed by Colgate and Johnson , 30 a significant portion of 
the difference between the observed and predicted light to medium nuclei ratio 
could be explained by a modulation mechanism which was rigidity dependent. 
However, the helium to medium ratio at low energies would then be very large 
and inconsistent with observations by large factors. Also, the energy density of 
the cosmic rays in the galaxy would then be very large. 

(G) The cosmic rays come from more than one source type. Assuming more 
than one source type in an equilibrium model does not resolve the difficulties 
unless one introduces a unique composition effect at the source. The case of a 
non-equilibrium model will be considered in the next section. 

IV. NON EQUILIBRIUM MODELS 

In exploring this possibility, it must be remembered that there is some 
evidence from meteorites that the cosmic rays have been constant on the aver- 
age within a factor of a few for a million years or more. The proposed sources 
must also produce the observed energy spectrum. The simplest, and in many 
respects the most likely, model of this type is one in which the particles from a 
single recent relatively close source are superimposed on the fluxes from all 
the other sources, whose effects can be represented reasonably well by the 
models discussed earlier. This possibility will be treated further in the next 
several paragraphs. 
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The potential path length distribution for particles observed from a single 
source will be very different from the ones already discussed. Assuming the 
particles are released within a short period, at least short compared to the time 
since their release , the path length of any observed particle will just be the velocity • 
time; hence, each energy interval will have a unique potential path length which 
will be proportional to the particle’s velocity. The relative intensity will be de- 
termined from an assumed diffusion picture, the work in sections HA and IIB, 
and the path length determined by the velocity, time, and interstellar density. 

In allowing this new degree of freedom, we shall return to what appears to 
uc tc be the most reasonable or at least the simplest model for the remainder 
of the cosmic rays, namely a uniform distribution of sources in space and 
time diffusing into an infinite isotopic medium. This model has at least two 
other advantages. First, the results are independent of the diffusion coefficient. 
Secondly, it applies both to the cosmic ray model in which the cosmic rays pervade 
the whole universe in equal intensity and the one in which the diffusion in the 
galaxy is slow, as present models would suggest, and there is a sufficient number 
of galactic sources spread fairly randomly — at least in the plane. The results 
of the earlier work showed that there was reasonable agreement with the experi- 
mental data at low energies, but not at high energies. 

If a relatively new close source has occurred, the high energy particles will 
be reaching us more efficiently than the lower energy ones even if the mean free 
path for scattering is independent of energy simply because the diffusion coeffi- 
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cient is proportional to ft. In order to make the discussion quantitative, isotropic 
diffusion will be assumed and specifically it will be assumed that the diffusion 
picture described by eq. (10) is valid. Eq. (10) gives the relative probabilities of 
various path lengths as a function of i and hence for ft in this case by using eq. 

(3) since (t - r ) is a constant now. The intensity of the various components for 
a given source spectral shape canthenbe obtained. There are effectively two 

- 3/2 

adjustable parameters in eq. (10), the constant multiplying ft at the front 
and the constant multiplying / 3 _1 in the exponential, since the other constants 
multiplying ft to give l can be absorbed in these two. The first of the constants 
was adjusted to give the observed light to medium nuclei ratio at high energies 
when this single source is added to the general background. At high energies 
this single source would represent about 40% of the total cosmic ray flux. The 
second constant was adjusted to give the best fit to the light to medium nuclei 
data. The choice of this second constant implies that the cosmic ray intensity 
will increase substantially in the future reaching a maximum of up to ten times 
the present value at high energies. Without knowing r 2 /% it is not possible to 
give aprecise time scale, but 10 5 to 10 6 years is probably reasonable. 

The results of the addition of this source calculated in the manner just 
described superimposed on a flux of isotropic sources in an isotropic medium 
is shown in figs. 12, 13, and 14. A reasonable agreement with experimental 

- 2.5 

data is possible with a source spectrum whose steepness lies between j s ~ W 
and j s ~ R~ 2 5 . 
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V. SUMMARY 

A comparison of the experimental data to the predictions of a large number 
of equilibrium models for the cosmic ray propagation showed that none of the 
models for which the mean free path was independent of rigidity yielded results 
which were in agreement with the experimental data, expecially the light to medium 
nuclei ratio. Further, the models selected were sufficiently varied and covered 
enough extremes so that other likely models which come to mind would give 
results intermediate between those predicted by some of the cases considered 
here, and, therefore, would be in disagreement with the experimental data also. 
as nas been nienll^eJ before, agree™"* +hp lip-ht to medium nuclei ratio 
can be obtained by choosing an appropriate dependence of the potential path length 
on energy. However, this dependence appears to be difficult to justify. 

One particular experimental result, namely the fluorine to oxygen ratio at 
low energies, about 100 MeV/nucleon, seems to be in basic disagreement with 
the fragmentation concept, and, therefore this cosmic ray abundance ratio and 
the relevant cross sections deserve serious attention. 

Several possible explanations for the difference between theoretical pre- 
dictions for the relative abundances of charges and the observations were 
considered with a particular emphasis on the light, medium, heavy and very 
heavy nuclei as groups, since considerable data exists for these species. It 
was seen that such considerations as the light nuclei being produced in a way 
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other than fragmentation, inadequate cross section data, different source spec- 
tra, a very complex galactic picture, poor experimental measurements, and 
solar modulation either did not seem to be likely explanations for the devia- 
tions on the basis of directly applicable data, or were very unlikely to provide a 
acceptable solution within the framework of our present understanding of other 
phenomena. 

One possible explanation which removes the difficulties associated with the 
light to medium nuclei ratio and still gives reasonable agreement with the medium, 
heavy, and very heavy relative abundances is a non-equilibrium model, wherein 
there has been one recent relatively close source that is making a significant 
contribution to the locally observed cosmic radiation. If this concept is correct, 
the composition of the high energy cosmic rays when examined in detail will re- 
flect the fact that about half of the high energy cosmic rays have gone through a rel- 
atively small amount of material. Further, it is possible that the energy spectrum of 
the "recent, close" source is slightly different than the average resulting in a 
variation in the composition at large cosmic ray energies. It would be very 
valuable, therefore, to have detailed composition measurements of the cosmic 
radiation at high energies, as well as low energies, to compare to the various 
possible predictions. 
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FIGURE CAPTIONS 


Fig. 1. The function 
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dx plotted as a function of — where — - 
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Fig. 2. The light to medium nuclei ratio for a potential path length of 4.5 g/cm 
for the source spectra indicated in the figure. 

Fig. 3. The helium to medium nuclei ratio for a potential path length of 4.5 
g/cm 2 for the source spectra indicated in the figure. 

Fig. 4. The heavy to medium and the very heavy to medium nuclei ratios for a 

2 2 5 

potential path length of 4.5 g/cm for source spectra of the form j s = R" 

Fig. 5. The light to medium nuclei ratio for the case of a point source ex- 
panding into an infinite isotropic medium for the source spectra indicated 
in the figure. 

Fig. 6. The helium to medium nuclei ratio for the case of a point source ex- 
panding into an infinite isotropic medium for the source spectra indicated 
in the figure. 

Fig. 7. The heavy to medium nuclei ratio and the very heavy to medium nuclei 
ratio for the case of a point source expanding into an infinite isotropic 

— 2 5 

medium for a spectrum of the form j = R 
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Fig. 8. The light to medium nuclei ratio for the case of a uniform distribution 
of sources over an infinite isotropic medium for the source spectra indicated 
in the figure. 

Fig. 9. The light to medium nuclei ratio for the case of a uniform distribution 
of sources contained within a sphere of radius R 0 for the source spectra 
indicated in the figure. 

Fig. 10. The helium to medium nuclei ratio for the case of a uniform distribu- 
tion of sources contained within a sphere of radius R 0 for the source 
spectra indicated in the figure. 

Fig. 11 . The light to medium nuclei ratio for a point source expanding into an 
infinite isotropic medium with K = a (1 + bR) for a source spectral of the 
form j g - R -2 ■ 5 . Note that fig. 5 shows that the choice of source spectra 
is not critical. 

Fig, 12. The light to medium nuclei ratio for the case of a relatively recent 
local source superimposed on the flux of isotropic sources in an isotropic 
medium described in the text in section III for the source spectral shapes 
given in the figure. References for data point are given in Reames and 
Fichtel 2 . 

Fig. 13. The helium to medium nuclei ratio for the case of a relatively recent 
local source superimposed on the flux to isotropic sources in an isotropic 
medium described in the test in section III for the source spectral shapes 
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given in the figure. References for data points are given in Fichtel and 
Reames. 2 

Fig. 14. The heavy and very heavy to medium nuclei ratios for the case of a 

relatively recent local source superimposed on the flux of isotropic sources 
in an isotropic medium described in the text in section III for the source 
spectral shapes given in the figure. References for data points are given 
in Fichtel and Reames, 2 except points indicated by squares, which are 
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from Omes and Webber. 
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